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Background. Middle East respiratory syndrome coronavirus (MERS-CoV) causes severe respiratory illness in humans. Funda¬ 
mental questions about circulating viruses and transmission routes remain. 

Methods. We assessed routinely collected epidemiologic data for MERS-CoV cases reported in Saudi Arabia during 1 January- 
30 June 2015 and conducted a more detailed investigation of cases reported during February 2015. Available respiratory specimens 
were obtained for sequencing. 

Results. During the study period, 216 MERS-CoV cases were reported. Full genome (n = 17) or spike gene sequences (n = 82) 
were obtained from 99 individuals. Most sequences (72 of 99 [73%]) formed a discrete, novel recombinant subclade (NRC-2015), 
which was detected in 6 regions and became predominant by June 2015. No clinical differences were noted between clades. Among 
87 cases reported during February 2015, 13 had no recognized risks for secondary acquisition; 12 of these 13 also denied camel 
contact. Most viruses (8 of 9) from these 13 individuals belonged to NRC-2015. 

Discussions. Our findings document the spread and eventual predominance of NRC-2015 in humans in Saudi Arabia during the 
first half of 2015. Our identification of cases without recognized risk factors but with similar virus sequences indicates the need for 
better understanding of risk factors for MERS-CoV transmission. 

Keywords. Middle East respiratory syndrome; MERS; coronavirus; MERS epidemiology; MERS transmission; MERS phylogeny; 
recombinant; Saudi Arabia. 


Middle East respiratory syndrome coronavirus (MERS-CoV) is 
known to cause severe respiratory illness in humans, with 
deaths recorded in 35%-40% of cases reported globally [1]. 
Since the first recognition of MERS in 2012, all cases reported 
to the World Health Organization have been linked to the Ara¬ 
bian Peninsula, with >85% of cases reported from Saudi Arabia 
[1]. Camels ( Camelus dromedarius) have been suspected as a 
reservoir for MERS-CoV, based on case investigations [2], sero¬ 
logic studies [3], and the isolation of virus from camels [4-8], 
Direct camel contact has also been identified as a risk factor 
for human illness [9]. Secondary human transmission has 
been demonstrated among close contacts of symptomatic 
cases, primarily following healthcare-associated exposures 
[10-12] and, to a lesser degree, household exposures [13]. 
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There is no definitive evidence of sustained human-to-human 
transmission in the community [14], 

MERS-CoV infection can exhibit a wide range of clinical 
manifestations, including mild or limited symptoms among 
those identified through contact tracing [11]. Prolonged viral 
shedding from the respiratory tract of those without obvious 
symptoms has been demonstrated [15], and transmission relat¬ 
ed to unrecognized cases has been suggested [12, 16] but not 
documented. 

MERS-CoV sequences obtained to date suggest periodic in¬ 
troductions of the virus into human populations, presumably 
from an animal reservoir, with subsequent limited chains of 
transmission in households and healthcare settings. The tempo¬ 
ral persistence of identified viral clades appears limited, consis¬ 
tent with an R 0 of <1 [17, 18]. Intervals between the beginning 
and end of the circulation of a clade vary, with longer intervals 
suggesting the existence of undetected human cases [19]. Cases 
and clusters continue to be reported from countries in or near 
the Arabian Peninsula, presenting an ongoing threat for broader 
transmission [20]. 

To assess the epidemiologic and clinical features of the dis¬ 
ease, we investigated all cases reported by the Saudi Arabia 
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Ministry of Health (MoH) during January-June 2015, and we 
attempted genetic sequencing on all available specimens. 

METHODS 

This investigation was part of an emergency public health re¬ 
sponse and was determined to be nonresearch by the MoH 
and Centers for Disease Control and Prevention (CDC) and 
therefore not subject to institutional review board review. 

Epidemiologic Investigation 

MERS-CoV Case Definition in Saudi Arabia 
At the time of this investigation, reporting in Saudi Arabia was 
required for all patients with clinical or radiologic evidence of 
MERS-CoV infection and a positive real-time reverse transcrip¬ 
tion-polymerase chain reaction (rRT-PCR) test result [21]. All 
rRT-PCR-positive cases identified at non-MoH facilities re¬ 
quired confirmation at MoH laboratories. 

January-June 2015 

We assessed the routinely collected epidemiologic information 
for all MERS-CoV cases reported by the MoH during 1 Janu- 
ary-30 June 2015, to provide a basic epidemiologic description. 
For this analysis, we included only individuals who met the case 
definition described above (ie, symptomatic cases). 

February 2015 

February 2015 was a period of increased reporting. To perform 
a more in-depth analysis, we collected additional information 
for all individuals with laboratory-confirmed MERS-CoV infec¬ 
tions during February 2015. This included all cases meeting the 
case definition as described above, as well as those identified as 
having a laboratory-confirmed case but no recognized symp¬ 
toms; individuals not meeting the case definition [21] were typ¬ 
ically identified through contact tracing. We reviewed available 
MoH case investigation records and data reported through the 
MoH Health Electronic Surveillance Network. We collected de¬ 
mographic information, medical history, outcome information, 
and treatment location. We assessed the likelihood of acquisi¬ 
tion from another person (secondary acquisition) by determin¬ 
ing whether a patient (1) was a healthcare professional (HCP), 
(2) had been admitted to a healthcare facility 2-14 days before 
illness onset, (3) had visited any healthcare facility in the 14 
days before illness onset, or (4) had direct contact with either 
another documented case of MERS-CoV infection or with 
someone with an acute respiratory illness of unknown cause 
in the 14 days prior to illness. When it was not possible to de¬ 
termine the criteria described above by using available informa¬ 
tion, we conducted telephone interviews (in Arabic) to collect 
additional exposure information. Proxies (a close friend or im¬ 
mediate family member who was familiar with the patient’s ac¬ 
tivities during this period) were interviewed if the case was 
deceased, still hospitalized, or too ill to participate. Among 
cases without any of the aforementioned risk factors for second¬ 
ary acquisition (hereafter referred to as sporadic cases), we 


asked during telephone interviews about the history of exposure 
to camels [9]. Interviewees were prompted to describe examples 
of camel exposures, including direct contact or visiting a live 
market, slaughterhouse, or race where camels were present. 
We also assessed travel history. 

Statistical Analysis 

Demographic and clinical characteristics were reported, and dif¬ 
ferences were assessed for significance by using % 2 , Wilcoxon 
rank sum, and Kruskal-Wallis tests, where appropriate. Data 
were analyzed using SAS, version 9.3 (SAS Institute, Cary, 
North Carolina). 

Molecular Detection and Sequencing 

Molecular testing was performed on all respiratory specimens 
available during January-June 2015. 

Specimens and Molecular Testing at the MoH 
Respiratory specimens, including nasopharyngeal and oral pha¬ 
ryngeal swabs, both separate and combined, nasopharyngeal and 
tracheal aspirates, and sputa collected from suspected MERS 
cases were tested at MOH laboratories by upE and ORFla 
rRT-PCR assays [22]. Available specimen aliquots (or RNA ex¬ 
tracts) that tested positive for MERS-CoV by both assays were 
shipped on dry ice to the CDC (Atlanta, Georgia) for sequencing. 

Molecular Testing at the CDC 

Sample aliquots (200-300 pL, if available) were extracted on a 
NucliSens EasyMAG (BioMerieux), and 100 pL of total nucleic 
acid elutes were recovered. The specimen extract were retested 
by MERS-CoV N2 and/or N3 rRT-PCR assays [23], and se¬ 
quencing was attempted on confirmed positive samples. Over¬ 
lapping nested primer sets were used for amplification and 
Sanger sequencing of the MERS-CoV spike genes and selected 
genomes (Supplementary Table 1). Amplicon sequencing was 
performed in both directions, using sequencing and internal 
amplification primers, with the BigDye Terminator v3.1 Cycle 
Sequencing Kit on a 3730x1 Genetic Analyzer (Thermo Fisher 
Scientific). Sequencher 5.3 software (Gene Codes) was used for 
sequence assembly and editing. 

Phylogenetic Analyses 

Nucleotide sequences were aligned using Clustal X, version 
1.83, implemented in BioEdit, version 7.2.5. Phylogenies were 
estimated using neighbor-joining and maximum likelihood 
methods implemented in Molecular Evolutionary Genetics 
Analysis, version 6.0621 [24], and Bayesian inference, using 
MrBayes v3.2.6 [25]. The neighbor-joining method used maxi¬ 
mum composite likelihood distance estimation and maximum 
likelihood used general time reversible (GTR) model of nucleo¬ 
tide substitution with y-distributed rate variation and a propor¬ 
tion of invariant sites (GTR + G + I). MrBayes was performed 
under a GTR model of nucleotide substitution with 4 categories 
of y-distributed rate heterogeneity and a proportion of invariant 
sites (GTR+ 4 +I). 
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NRC-2015 


-KT806020 1 |Hu/Riyadh-KSA-10182/201515-2015 • 

- KT225476.2|MERS-CoV/THA/CU/17_06_2015|6-2015 

- KT006149.1 |Hu/ChinaGD0111 -2015 

KT806039.1 |Hu/Hufuf-KSA-12354/2015|6-2015 • 

—KT806037.1 |Hu/H ufuf-KSA-12177/2015|6-2015® 

—KT806032.1 |Hu/Hufuf-KSA-11765/2015|5-2015 • 
KT806041 11Hu/Hufuf-KSA-11078/2015|5-2015 • 

—KT806040.1 |Hu/Hufuf-KSA-12597/2015|6-2015 • 

KT806036.11Hu/Hufuf-KSA-12124/2015|6-2015 • 

KT806035.1 |Hu/Hufuf-KSA-12055/2015|6-2015 • 

KT806034.1 |Hu/Hufuf-KSA-11876/2015|5-2015 • 

KT806031 1 |Hu/Hufuf-KSA-11767/2015|5-2015 • 
KT806030.1|Hu/Hufuf-KSA-11853/2015|5-2015# 
KT806029.1 IHu/Hufuf-KSA-11869/2015|5-2015 • 

KT806028 1|Hu/Hufuf-KSA-11905/2015|5-2015 • 
KT806026.1 |Hu/Hufuf-KSA-11476/2015|5-2015 • 
KT806025.1 |Hu/Hufuf-KSA-11442/2015|5-2015 • 

KT806024.1 |Hu/Hufuf-KSA-11401/2015|5-2015 • 

KT806033.1 |Hu/Hufuf-KSA-10357/2015|4-2015 • 
KT805990.1 |Hu/Hufuf-KSA-11003/2015|5-2015 • 

KT806046 1 |Hu/Hufuf-KSA-11002/2015|5-2015 • 
KT805989.1 |Hu/Hufuf-KSA-11150/2015|5-2015 • 

KT805987.1 |Hu/Hufuf-KSA-11080/2015|5-2015 • 

—KT806027.1 |Hu/Hufuf-KSA-11882/2015|5-2015# 
KT806038.1 |Hu/Hufuf-KSA-12225/2015|6-2015 • 

KT029139.1 |MERS-CoV/KOR/KNIH/002_05_201515-2015 
- KT806049.1 |Hu/Riyadh-KSA-3181 /2015|2-2015 • 

-KT806016.1 |Hu/Duba-KSA-313/2015|4-2015# 

KT805996.1 |Hu/Najran-KSA-C22374/2015|3-2015 • 

KT806015.1 |Hu/Makkah-KSA-907/2015|4-2015 • 

—KT806054 1 |Hu/Najran-KSA-C20915/2015)2-2015 • 
iKT806048 1 |Hu/Khobar-KSA-6736/2015|2-2015 • 

KT805984.1 |Hu/Dhahran-KSA-6740/2015|2-2015 • 

KT805993.1 |Hu/Jeddah-KSA-C21760/2015|3-2015# 
KT806045.1|Hu/Jeddah-KSA-C21271/2015|2-201 5# 

KT806044.1 |Hu/Jeddah-KSA-C20843/2015|2-201 5# 

KT805982.1 |Hu/Hufuf-KSA-6648/2015| 1 -2015# 

KT806053.1 |Hu/Kharj-KSA-2598/2015|2-201 5# 

— KT0264 54.1 |Hu/Riyadh_KSA_4050_2015|3-2015 • 

— KT806047.1|Hu/Hufuf-KSA-9158/2015|3-201 5# 

KT806052.1 |Hu/Kharj-KSA-2599/2015|2-2015 • 

KT806050.1 |Hu/Riyadh-KSA-3065/2015|2-2015 • 
iKT805978.1 |Hu/Riyadh-KSA-3550/2015|2-2015 • 

1 KT806003.1 |Hu/Rass-KSA-3318/2015|2-2015 • 

—KT805983.1 |Hu/Hufuf-KSA-6674/2015|2-2015 • 

H KT806014.1 |Hu/Riyadh-KSA-7867/2015)4-2015 • 

KT806023.1 |Hu/Aflaj-KSA-10802/2015|5-2015 • 

KT806002.1 |Hu/Buraydah-KSA-3317/2015|2-2015 • 

KT806012.1 |Hu/Riyadh-KS A-7591/2015|3-2015 • 

KT806022.1|Hu/Ounfidhah-KSA-Rs1338/2015|5-2015# 

KT805967.1 1 Hu/Riyadh-KS A-3331 /2015|2-2015# 

_M KT805985 1 |Hu/Hufuf-KSA-8271 /2015|3-2015 • 

_| ' KT805986.1 |Hu/Hufuf-KSA-8527/2015)3-2015 • 

■KT805968.1 |Hu/Riyadh-KSA-3160/2015|2-2015 • 

KT806009.1 |Hu/Riyadh-KSA-5390/2015|3-2015 • 

KT805964 1 |Hu/Riyadh-KSA-2600/2015|2-2015 • 

KT805991 1 |Hu/Najran-KSA-Rs1207/2015)5-2015# 

KT806018.1 1 Hu/Riyadh-KS A-10131 /2015|5-2015 • 

KT805980.1 |Hu/Riyadh-KS A-3487/2015)2-2015 • 

KT026453.1 |Hu/Riyadh_KSA_2959_2015|2-2015 • 

KT805979.1 |Hu/Riyadh-KSA-3486/2015|2-2015 • 

|—KT806008.1 |Hu/Shaqra-KSA-4699/2015)3-2015# 

KT806017.1 |Hu/Taif-KSA-15167/2015|5-2015 • 

KT806007.1 |Hu/Riyadh-KSA-4606/2015|3-2015 • 

_ KT806011.1 |Hu/Riyadh-KSA-7244/2015)3-2015 • 

KT806005.1 |Hu/Riyadh-KSA-3782/2015|2-2015 • 

KT805974.1 |Hu/Riyadh-KSA-4071 /2015|2-2015 • 

KT805965.1 |Hu/Riyadh-KS A-2622A/2015)2-2015 • 

KT806051.1 |Hu/Riyadh-KSA-2716/2015|2-2015 • 

—KT806021.1 |Hu/Riyadh-KSA-6966/2015|3-2015# 

KT805977.1 1 Hu/Riyadh-KS A-3562/2015|2-2015 • 

KT805971.1 |Hu/Riyadh-KSA-3123/2015|2-2015 • 

_ KT805976.1 |Hu/Riyadh-KSA-3929/2015|2-2015• 

—KT805975.1 |Hu/Riyadh-KSA-4067/2015|2-2015 • 

—KT806010.1 |Hu/Riyadh-KSA-7179/2015|3-2015 • 


Figure 1 . Phylogeny of Middle East respiratory syndrome coronavirus (MERS-CoV) spike sequences. Phylogenetic analysis of the MERS-CoV spike gene coding region 
(4062 nucleotides) obtained from 99 cases in this study and 251 previously published sequences available in GenBank. Major clades A and B are indicated by vertical bars. 
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Genetic Recombination and Ancestral Analysis 

Putative recombination events were identified using Recombi¬ 
nation Detection Program software, version 4.70 (RDP4; avail¬ 
able at: http://web.cbio.uct.ac.za/~darren/rdp.html), with the 
default settings [26]. The complete genome sequence of each of 
the viruses in the NRC-2015 clade was aligned with the ge¬ 
nomes outside the clade. The multiple sequence alignment was 
then imported into the RDP software for detection of recombi¬ 
nation. The software uses several algorithms, including GENE- 
CONV, BootScan, MaxChi, Chimaera, SiScan, and 3Seq, to 
detect putative recombination events. The potential minor 
and major parental sequences and the beginning and end 
breakpoints of the potential recombinant sequences were also 
defined by RDP4 software. Putative recombinant events were 
considered significant when a P value of < .05 was observed 
for the same event, using >4 algorithms. 

Time estimates to the most recent ancestor were calculated 
using the Bayesian Markov Chain Monte Carlo (MCMC) meth¬ 
od implemented in BEAST vl.8.2 [27], The coding regions 
(ORFlab, S, ORF3, ORF4a, ORF4b, ORF5, E, M, and N) in 
the genomes grouping within NRC-2015 were concatenated, 
and the E1KY+ T 4 substitution model was used with indepen¬ 
dent rates for each of the positions in the codon. A lognormal 
relaxed molecular clock (uncorrelated) was used with Gaussian 
Markov random field Bayesian skyride coalescent. Bayesian 
MCMC analysis was run for 25 million steps. Parameters for 
tMRCA, rate, and trees were sampled every 5000 steps, with 
the first 10% removed as burn-in. Time estimate values thus ob¬ 
tained were also compared with strict and exponential relaxed 
clock models. 

RESULTS 

During 1 January-30 June 2015, 216 MERS-CoV cases from 10 
of the 13 regions of Saudi Arabia were reported by the MoH; 
MERS-CoV-positive individuals with no recognized symp¬ 
toms, and who therefore did not meet the case definition, were 
not included. The longest period between case reports was 11 
days. Among these 216 cases, 214 were hospitalized, and 102 
(47%) died. Most patients were male (161 [75%]) and of Saudi 
nationality (147 [68%]). Median age was 56 years (range, 20-93 
years). 

Molecular Analysis of MERS-CoV Strains 

Of the 216 symptomatic cases reported during the study period, 
124 had respiratory specimens available for further testing at the 
CDC; 1 specimen was also available from an individual with no 
recognized symptoms who did not meet the case definition. Of 


the 125 available respiratory specimens collected during 6 
January-3 June 2015, spike gene sequences were obtained 
from 99 (Supplementary Table 2). Phylogenetic analysis of 
the 99 spike sequences with 251 previously published sequenc¬ 
es (Figure 1) revealed a novel subclade (henceforth referred to 
as NRC-2015) within clade B [18]; this subclade contained 
most (72 of 99 [72.7%]) of the 2015 sequences and included 
sequences recently reported from travel-associated MERS cases 
in the Republic of Korea (accession number KT029139), China 
(KT006149), and Thailand (KT225476). The neighbor joining 
tree bootstrap value supporting NRC-2015 was low (<70%), 
which may reflect the recombinant character of the clade se¬ 
quences (see below) or the lower number of informative sites 
in the spike gene. Elowever, similar tree topologies reproducing 
NRC-2015 were obtained with the spike gene sequences by 
both maximum likelihood and Bayesian analyses, supporting 
this finding (data not shown). To more fully characterize NRC- 
2015, 17 MERS-CoV-positive samples with sufficient available 
volume were selected for genome sequencing, including 12 
NRC-2015 sequences from different time points and 5 from out¬ 
side of the subclade. 

Phylogenetic analysis of these full genomes with genome se¬ 
quences available on GenBank reproduced NRC-2015 with 
high support values (Supplementary Figure 1). Mutational anal¬ 
ysis of the 12 new NRC-2015 genomes with previously pub¬ 
lished NRC-2015 genomes obtained from human samples 
collected in Saudi Arabia (n = 2), Korea/China (n = 2), and 
Thailand (n = 1) revealed 9 unique nucleotide substitutions 
common to all NRC-2015 members, of which 3 were nonsy- 
nonymous, with predicted amino acid changes in the ORFla 
(n = 2) and matrix (n = 1) proteins (Supplementary Table 3). 
Recombination analysis on the newly available genome se¬ 
quences from NRC-2015 identified 2 possible recombination 
events involving sequences from outside the clade as potential 
minor and major parental strains. The first event had a predict¬ 
ed breakpoint at nucleotide position 17 475 (99% confidence in¬ 
terval [Cl], 13 502-19 074), located in ORFlab, and the second 
event had a predicted breakpoint at 23 976 (99% Cl, 23 571- 
24 862), located in the spike gene. Recombination analysis 
was performed using RDP software, and events detected with 
a P value of < .05 were considered evidence of true recombina¬ 
tion (Supplementary Table 4). 

To date the emergence of NRC-2015, MCMC analysis was 
performed on the concatenated coding regions of the genomes 
grouping within NRC-2015, using BEAST. The most recent 
common ancestor of the virus was approximately 0.85 years 


Figure 1 continued. The novel subclade (NRC-2015) is repositioned for clarity. Spike gene sequences obtained in this study are marked with solid circles; red circles 
indicate samples with available genome sequences (Supplementary Figure 1 ). The tree was constructed by the neighbor-joining method, and bootstrap resampling values 
(1000 replicates) >70% are indicated above the respective nodes. Numbers in brackets following some strain identifiers are the number of identical sequences with the 
same location and sample collection time. The scale bar shows the genetic distance as nucleotide substitutions per site. To compare persistence and geographic dis¬ 
tribution of NRC-2015 with other viruses identified during 2012-2015, we assigned virus sequences to different subclades (depicted in color), based on previously de¬ 
scribed clades [19, 28-30], 
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Figure 2. Middle East respiratory syndrome coronavirus (MERS-CoV) cases reported by the Saudi Arabia Ministry of Health (MoH) during 1 January-30 June 2015. A, MERS- 
CoV cases, by clade and week of illness onset. B, Sequenced viruses, by clade and month of illness onset. The 216 cases reported by the MoH during 1 January-30 June 2015 
are included. These graphs do not include the single sequence from the individual with no recognized symptoms who did not meet the case definition. NRC-2015 was defined 
using the spike gene phylogenies (Figure 1). 


(95% Cl, .45-1.34) before sample KT225476.2 (17 June 2015), 
or August 2014 (Supplementary Table 5). 

Epidemiologic Investigation 

Among the 216 cases reported during 1 January-30 June 
2015, NRC-2015 was first detected in a case with onset in mid- 
January 2015 (Figure 2A). During the study period, NRC-2015 
viruses were detected in 6 regions of Saudi Arabia (Figure 3), 
and the proportion of patients identified with NRC-2015 in¬ 
creased steadily over time (Figure 2 B). 

NRC-2015 was next compared to past and present subclades 
within clade B, using sequences available in GenBank (Figure 4). 
NRC-2015 was more widely distributed geographically than any 


other identified members of clade B. The duration of circulation 
of recognized subclades ranged from 16 to 665 days. At the con¬ 
clusion of our investigation period, NRC-2015 had been circu¬ 
lating for 135 days, which was longer than 7 of 9 other identified 
subclades. In our analysis, the longest circulating subclade re¬ 
ported was Riyadh_KKUH-l_2014, which was first detected 
in July 2013 and was still circulating as of May 2015. During 
our investigation period, 10 of 99 sequenced viruses belonged 
to Riyadh_KKUH-l_2014. No viruses belonging to clade A 
were detected. 

A comparison of patients infected with NRC-2015 versus 
other circulating viruses revealed no significant differences in 
age, sex, rate of mortality, time between onset of symptoms 
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Figure 3. Map of Middle East respiratory syndrome coronavirus (MERS-CoV) cases, by subclade and region of Saudi Arabia. The 216 cases reported by the Saudi Arabia 
Ministry of Health during 1 January-30 June 2015, are included. This map does not include the single sequence from the individual with no recognized symptoms who did not 
meet the case definition. NRC-2015 was defined using the spike gene phylogenies (Figure 1). 


and death, or length of hospital stay among fatal cases (Supple¬ 
mentary Table 6). There was also no difference in mean cycle 
threshold values, a proxy for virus load, with respiratory speci¬ 
mens containing NRC-2015 versus other clades, although these 
were not adjusted for timing of specimen collection (Supple¬ 
mentary Table 6). 

For our more detailed analysis of cases reported during 1-28 
February 2015, we identified 87 MERS-CoV-positive patients 
(Table 1). Of these, 77 patients (89%) satisfied the case defini¬ 
tion for routine reporting and required hospitalization; the re¬ 
maining 10 individuals (11%) had no recognized symptoms 
(and did not satisfy the case definition) but are included in 
this analysis. The 87 patients were reported from 35 different 
healthcare facilities across 7 regions in Saudi Arabia; 17 of 
these facilities reported >2 cases within the same 14-day period. 
Of these 87 patients, sequences could be obtained from 34, of 
which 24 (71%) were associated with NRC-2015. No clinical 
differences were apparent when comparing NRC-2015 to 
other circulating viruses (Table 1). 

The 87 patients with laboratory-confirmed disease reported 
during February were also classified according to their reported 
exposures during the 2 weeks before illness onset. Record review 
and interviews were conducted during 11-25 March 2015. 


Among the 87 cases, 51 were classifiable using information ob¬ 
tained by the initial case investigation. Interviews were attempt¬ 
ed for the remaining 36 patients. Of these, 28 (78%) were 
interviewed; 1 individual refused to participate, and 7 patients 
were not available. Proxy interviews were conducted for 22 of 
28 interviews, including for 18 patients who were deceased 
and for 4 of 10 patients who survived. 

Among the 87 patients, 13 (15%) were determined to have 
had household contact with a confirmed MERS-CoV case, 14 
(16%) were FICPs, 21 (24%) were inpatients in a healthcare fa¬ 
cility, 16 (18%) were hospital visitors, and 10 (11%) were unable 
to be classified owing to a lack of available information (Table 1). 
Notably, 13 patients (15%) denied exposure to a healthcare fa¬ 
cility or to a person with acute respiratory illness in the 2 weeks 
before illness onset and were classified as sporadic cases (Ta¬ 
bles 1 and 2); among these, 1 individual reported visiting a 
camel farm in the 2 weeks before illness onset. Among the 13 
sporadic cases, 2 were available for interview, and 11 were inter¬ 
viewed by proxy. Among the 11 interviewed by proxy, 9 were 
deceased and 2 were too ill to participate in the interview. Se¬ 
quences were obtained for 9 sporadic cases, and 8 (89%) were 
NRC-2015, including the individual who had visited the camel 
farm. 
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A 


• Riyadh_3 

• Hafr-AI-Batin_l 

• Riyadh KKUH-1 2014 

• Riyadh KKUH-2 2014 

• NRC 2015 


• Al-Hasa 

• Buraidah_l 

• Jeddah KFH 

• Riyadh PSMMC 2015 

• Jeddah-1 2015 


• • • 


• # •• # 

1 Jan 2013 13 Apr 2013 24Jul2013 3 Nov 2013 13 Feb 2014 26 May 2014 5 Sep 2014 16 Dec 2014 28 Mar 2015 8Jul2015 

Specimen Collection Date 




Summary of clade dates and locations 


Clade 

Date Range 

Day, No. 

Distinct 
Locations, No. 

Riyadh_3 

Feb 2013-Nov 2013 

275 

4 

Al-Hasa 

Apr 2013-Jun 2013 

63 

1 

Hafr-AI-BatinJ 

Jun 2013-Sep 2013 

100 

3 

Buraidah_1 

Jun 2013-Aug 2013 

55 

2 

Riyadh KKUH-1 2014 

Jul 2013-May 2015 

665 

7 

Jeddah KFH 

Mar 2014-Apr 2014 

33 

3 

Riyadh KKUH-2 2014 

Mar 2014-Apr 2014 

33 

1 

Riyadh PSMMC 2015 

Jan 2015-Jan 2015 

21 

1 

NRC-2015 

Jan 2015-Jun 2015 

135 

15 

Jeddah-1 2015 

Mar 2015-Apr 2015 

16 

2 


Figure 4. Geotemporal distribution of Middle East respiratory syndrome coronavirus (MERS-CoV) clades in Saudi Arabia during 1 January 2013-30 June 2015. For human- 
derived sequences in each subclade within the broader clade B (as defined using the spike gene phylogenetic analysis in Figure 1), we used GenBank to determine length of 
persistence, using the earliest and most recent date of detection (A), and we determined geographic distribution, based on city of detection within Saudi Arabia (B). a 2015 
depicts clades detected through 30 June 2015. 


DISCUSSION 

Our study provides the first nationwide epidemiologic and phy¬ 
logenetic description of MERS-CoV in humans in Saudi Arabia, 
allowing for a better understanding of transmission during a 


period of continuous, frequent reporting of new cases. Our in¬ 
vestigation demonstrates the emergence, persistent circulation, 
and eventual predominance of a recently identified recombi¬ 
nant MERS-CoV clade in humans that includes viruses detected 
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Table 1. Epidemiologic and Clinical Characteristics of Patients Identified 
With Laboratory-Confirmed Middle East Respiratory Syndrome Coronavirus 
Infection in February 2015 in Saudi Arabia 


Characteristic 

Sequenced 

NRC-2015 Other 

Clade, Clades, 

No. (%) No. (%) 

(n = 24) (n = 10) 

No 

Sequence, 
No. (%) 

(n = 53) 

Total, 
No. (%) 
(n = 87) 

Sex 

Male 

17 (71) 

6 (60) 

39 (74) 

62 (71) 

Female 

7(29) 

4(40) 

14(26) 

25 (29) 

Age group, y 

18-34 

2 (8) 

1 (10) 

13 (25) 

16 (18) 

35-64 

14 (58) 

9 (90) 

24 (45) 

47 (54) 

>65 

8 (33) 

0(0) 

15 (28) 

23 (26) 

Unknown 

0(0) 

0(0) 

1 (2) 

1 (1) 

Nationality 

Saudi 

14 (58) 

6 (60) 

36 (68) 

56 (64) 

Non-Saudi 

10 (42) 

4(40) 

17 (32) 

31 (36) 

Underlying illness 

Any 

14(58) 

7(70) 

30 (57) 

51 (59) 

Diabetes 

10 (42) 

4(40) 

16 (30) 

30 (34) 

Heart disease 

1 (4) 

1 (10) 

3(6) 

5(6) 

Hypertension 

10 (42) 

1 (10) 

16 (30) 

27 (31) 

Exposure classification 

Household 

1 (4) 

2 (20) 

10 (19) 

13 (15) 

Healthcare provider 

4(17) 

1 (10) 

9(17) 

14 (16) 

Inpatient 

3 (13) 

3 (30) 

15 (28) 

21 (24) 

Hospital visitor 

6 (25) 

2 (20) 

8(15) 

16 (18) 

Sporadic 8 

8(33) 

1 (10) 

4(8) 

13 (15) 

Undetermined 13 

2 (8) 

1 (10) 

7(13) 

10 (11) 

Hospitalized 

No 

0 (0) 

0(0) 

10 (19) 

10(11) 

Yes, survived 

10 (42) 

4(40) 

22 (42) 

36 (41) 

Yes, died 

14 (58) 

6(60) 

21 (40) 

41 (47) 


a Patients who were not a contact of a case or a person with acute respiratory illness and did 
not report exposure to a healthcare facility in the 2 weeks prior to illness onset. 
b Patients for whom sufficient data were not available. 


in the Republic of Korea [31, 32], Thailand (accession number 
KT225476), and China in 2015 [31, 33]. Previous documenta¬ 
tion of the duration of circulation in humans of 4 different 
MERS-CoV clades in Saudi Arabia during 2012-2013 noted 
an average detection time of 98 days [19], In contrast, we dem¬ 
onstrate that NRC-2015 has persisted longer than most previ¬ 
ously documented clades. NRC-2015 was found to eventually 
predominate over the 6-month study period and attain a wide 
geographic distribution in a comparatively short period. While 
this apparent emergence and clade displacement is suggestive of 
greater epidemiologic fitness [34], we observed no clinical dif¬ 
ferences between NRC-2015 and other clades; the implications 
for virus replication and transmission need further study. Dur¬ 
ing preparation of this manuscript, sequences obtained from 
camels in Oman in May 2015 [35] and Saudi Arabia during 
July 2014-April 2015 [36] were reported that showed similar 


recombination features and phylogenetic association with 
NRC-2015. In camels, NRC-2015 (referred to as lineage 5 
[36]) was first detected in July 2014 and became predominant 
in Saudi Arabia during a period that overlaps with our study, 
corroborating our findings of an increased prevalence in hu¬ 
mans relative to other clades. 

Recombination has been documented among Co Vs [37] and 
has been linked to the emergence of more-pathogenic strains of 
some animal CoVs [38-40]. Evidence of intraspecies recombi¬ 
nation has also been found with the human Co Vs HKU1 [41], 
NL63 [42],OC43 [43], and, more recently, MERS-CoV [44]. Ge¬ 
nome analysis of human MERS-CoV strains from Saudi Arabia 
in 2015 and the recent outbreak in South Korea/China [31-33] 
and camels as noted above [35, 36] revealed a probable signature 
recombination event between 2 different parental clade B virus¬ 
es involving a region of the ORFlab and spike genes. We con¬ 
firmed this finding and documented an increasing prevalence of 
this virus in humans among samples collected since January 
2015 from geographically distant communities in Saudi Arabia. 
Similar to recent reports [33], we estimate that this recombinant 
virus emerged sometime in mid-to-late 2014. Based on recently 
available sequence data from camels in Saudi Arabia, NRC-2015 
(lineage 5) was predicted to have diverged between December 
2013 and June 2014 [36], 

In our study, further analysis of 87 patients with laboratory- 
confirmed MERS-CoV reported in February revealed 13 indi¬ 
viduals with no recognized risks for secondary acquisition; 
none of these 13 reported direct camel contact, although 1 in¬ 
dividual reported visiting a camel farm. Of those sequenced, 
most were infected with genetically very similar viruses, sug¬ 
gesting a potential for limited transmission from those with un¬ 
recognized MERS-CoV infection. These findings highlight the 
importance of strengthened epidemiologic and laboratory 
surveillance. 

Most cases identified in Saudi Arabia in February had docu¬ 
mented exposure to healthcare facilities, a well-demonstrated 
risk factor for MERS-CoV infection [10-12]. Seventeen of 35 
affected facilities in Saudi Arabia in February experienced 
MERS-CoV infection clusters. Moreover, 16 of 87 patients in 
February (18%) were visitors to healthcare facilities. This is sim¬ 
ilar to the 2014 Jeddah outbreak, where 17% of investigated 
cases were visitors [11]. Recommendations to limit visitation 
in facilities with ongoing MERS-CoV transmission should be 
reinforced to limit these exposures. 

Our investigation, which was performed as part of an emer¬ 
gency public health response, is subject to several limitations. 
First, specimens were not available for all cases during the 
study period, meaning that many viruses remained untyped; 
however, we observed no demographic differences between 
cases who had specimens sequenced and those who did not. 
Second, since its emergence in 2012, surveillance and sequenc¬ 
ing of MERS-CoV strains has been incomplete; variations in 
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Table 2. Sporadic Middle East Respiratory Syndrome (MERS) Coronavirus Cases Reported in Saudi Arabia During February, 2015 


Patient 

Interview 

Clade (GenBank 3 ) 

Region 

Age, y 

Sex 

Saudi 

Nationality 

Survived 

Underlying 

Illness 

Ever 

Smoker 

Occupation 

1 

Proxy 

NRC-201 5 g (KT806053) 

Riyadh 

56 

F 

Yes 

No 

Yes 

Unknown 

Unemployed 

2 

Self 

NRC-201 5 g (KT806052) 

Riyadh 

26 

M 

Yes 

Yes 

Yes 

Yes 

Security 

3 

Proxy 

NRC-2015 G (KT806048) 

Eastern 

49 

M 

Yes 

No 

Yes 

Yes 

Sales 

4 

Proxy 

NRC-201 5 g (KT806044) 

Jeddah 

73 

F 

No 

No 

Yes 

No 

Unemployed 

5 

Proxy 

NS 

Madinah 

70 

M 

Yes 

No 

Yes 

No 

Retired 

6 

Proxy 

NS 

Riyadh 

79 

M 

Yes 

No 

Yes 

Yes 

Unemployed 

7 

Proxy 

NRC-2015 s (KT806050) 

Riyadh 

49 

F 

Yes 

Yes 

Yes 

No 

Unemployed 

8 

Proxy 

NRC-201 5 g (KT806054) 

Najran 

47 

M 

No 

No 

Yes 

No 

Blacksmith 

9 

Proxy 

NS 

Eastern 

41 

M 

Yes 

Yes 

No 

No 

Unemployed 

10 

Proxy 

Other s (KT805966) 

Riyadh 

51 

M 

Yes 

No 

Yes 

Yes 

Retired 

11 

Proxy 

NRC-201 5 g (KT806045) 

Jeddah 

65 

F 

No 

No 

Yes 

No 

Unemployed 

12 

Proxy 

NS 

Qasim 

45 

M 

No 

No 

Yes 

No 

Joiner/blacksmith 

13 b 

Self 

NRC-2015 s (KT806002) 

Qasim 

42 

M 

No 

Yes 

No 

Yes 

Driver 


In the 2 weeks before becoming ill, none had known contact with a patient with MERS, a person with severe respiratory illness, and anyone mildly ill, and none reported working in or visiting a 
healthcare facility. 

Abbreviations: G, full genome sequence obtained; S, spike region sequence obtained; NS, not sequenced. 
a GenBank accession number. 
b Visited a camel farm. 


sequence availability and documentation might have influenced 
the extent of persistence and geographic spread that we have 
determined for past circulating virus strains. Case definitions, 
testing practices, and testing locations have also changed during 
this period. Third, although we were able to obtain full genome 
sequences from 12 NRC-2015 samples, all of which possessed the 
expected recombinant signal, our sequencing was mostly limited 
to the spike gene alone, which poses the risk of misclassifying re¬ 
combinant viruses as belonging to NRC-2015. This is illustrated 
in the recent study by Sabir et al [36], which reported multiple 
novel recombinant viruses in camels, including recombinants 
between NRC-2015 (lineage 5) and other virus clades. Fourth, be¬ 
cause of the high morbidity and mortality of MERS-CoV infection, 
interviews with cases were not always possible, necessitating the 
use of proxies. It is possible that, combined with issues of recall, 
the quality of the information collected varied. Of particular 
consideration, 11 of 13 sporadic cases were classified on the 
basis of interviews with proxies, and pre-illness exposures 
might not have been accurately recognized and reported. Fifth, 
some camel exposures may have gone unrecognized because of 
disincentives for reporting camel exposures, given their cultural 
and economic significance in Saudi Arabia. Sixth, given the ex¬ 
isting evidence of association between MERS-CoV illness and 
pre-illness healthcare exposure or exposure to sick individuals 
[10, 11, 13], our risk classification was hierarchical; that is, re¬ 
ported exposure to a setting where secondary acquisition was 
likely took precedence over reported exposure to camels. As 
such, we did not assess camel exposures in individuals with 
recognized risks for secondary acquisition. 

Finally, although we have attempted to link the results of 
our epidemiologic investigation with MERS-CoV sequences 


obtained from investigated cases, we cannot fully assess the pos¬ 
sible role of virus introductions from nonhuman sources. Re¬ 
cent phylogeny of MERS-CoV sequences from camels in 
Saudi Arabia indicated that the novel recombinant subclade (re¬ 
ferred to as NRC-2015 in our manuscript) was also predomi¬ 
nant in camels during a period overlapping with our study 
[36]. As such, our detection of closely related viruses in humans 
might in part reflect multiple introductions from camels with 
similar strains. Virus introductions from other currently un¬ 
identified sources might also be factor. Virus transmission dy¬ 
namics within and between human and nonhuman sources of 
MERS-CoV will likely influence transmission routes in ways not 
yet fully understood. 

This investigation describes the emergence, persistence, and 
widespread circulation of a novel recombinant MERS-CoV in 
Saudi Arabia. A lack of clearly defined epidemiologic links in 
some cases highlights the need for ongoing intensive epidemiolog¬ 
ic and laboratory surveillance to better understand MERS-CoV 
transmission and to focus infection prevention and control efforts. 
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